Mapping the Berry Curvature from Semiclassical Dynamics in Optical Lattices 
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We propose a general method by which experiments on ultracold gases can be used to determine 
the topological properties of the energy bands of optical lattices, as represented by the map of the 
Berry curvature across the Brillouin zone. The Berry curvature modifies the semiclassical dynamics 
and hence the trajectory of a wave packet undergoing Bloch oscillations. However, in two dimensions 
these trajectories may be complicated Lissajous-like figures, making it difficult to extract the effects 
of Berry curvature in general. We propose how this can be done using a "time-reversal" protocol. 
This compares the velocity of a wave packet under positive and negative external force, and allows a 
clean measurement of the Berry curvature over the Brillouin zone. We discuss how this protocol may 
be implemented and explore the semiclassical dynamics for three specific systems: the asymmetric 
hexagonal lattice, and two "optical flux" lattices in which the Chern number is nonzero. Finally, we 
discuss general experimental considerations for observing Berry curvature effects in ultracold gases. 

PACS numbers: 03.65.Sq, 03.65. Vf, 67.85.-d 
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I. INTRODUCTION 

One of the most interesting and surprising develop- 
ments in the band theory of solids was the realization 
that the physical properties of any resulting energy band 
are not only encoded in its energy spectrum, e(k), for 
all wavevectors k in the Brillouin zone (the "band struc- 
ture" in the usual sense). In addition, there are impor- 
tant physical consequences related to the topology of the 
eigenstates that form the band yj, |2( . 

The importance of topological features of the energy 
eigenstates was first pointed out in the seminal work of 
Thouless et al. Q in the context of the integer quantum 
Hall effect. This work showed that, for two-dimensional 
(2D) lattices, the set of energy eigenstates that form 
an isolated band are characterized by an integer-valued 
Chern number, C, which can be nonzero when time- 
reversal symmetry is broken. The Chern number is a 
topological invariant of the band: its value cannot be 
changed by continuous evolution of the physical system 
without closing the gap to another band. Furthermore 
Ref. showed that the Chern number has direct physical 
consequences: a filled band has a quantized Hall con- 
ductance equal to C times the quantum of conductance, 
corresponding to the existence of C chiral edge modes. 

Soon after the work of Thouless et al.\^, the concept 
of the Berry phase was formulated The integer in- 
variant of Ref. was quickly interpreted @ in terms of 
the integral of the Berry curvature in momentum space, 
f2(k), over the Brillouin zone. This Berry curvature, 
f2(k) (defined below), is a gauge-invariant property of the 
band structure which is predicted to have direct physical 
consequences when the band is partially filled [a, 0] • A 
full description of the properties of a nondegenerate en- 
ergy band therefore requires a specification of both the 
spectrum, e(k), and the Berry curvature, Q(k). For sys- 
tems with additional global symmetries (time-reversal, 
particle-hole, or chiral) other forms of topological invari- 
ant can appear [l], 0. 

These considerations have become of immediate im- 



portance in the field of ultracold gases in optical lat- 
tices. While the Berry curvature vanishes for the simplest 
forms of optical lattice, recent experiments have stud- 
ied more complicated 2D optical lattices |8hT(| for which 
simple variants exist in which the (local) Berry curvature 
is nonzero. Furthermore, there exist several theoretical 
proposals for forms of optical lattice in which the neu- 
tral atoms feel an effective magnetic field with high flux 
density [TM3- 

These proposals are of great interest, 
as they offer the opportunity to study ultracold atoms 
in strong effective magnetic fields. These lattices break 
time-reversal symmetry in a way that leads to bands with 
nonzero Chern numbers [lj, lll5Hl7 1 • A key motivation of 
this work is to propose how the bandstructures of these 
lattices might be characterized experimentally. We shall 
focus on the "optical flux lattices" proposed in Refs. [l6l 
and [III 

A particle picks up a Berry phase when it travels 
adiabatically around a closed contour. Here, we are 
concerned with contours in two-dimensional momentum 
space, with k in the Brillouin zone. As usual, the Berry 
phase can be expressed as the integral of the Berry cur- 
vature, f2(k), over the surface bounded by the contour. If 
the integral is over the entire Brillouin zone, the periodic- 
ity in k space means that the closed contour is equivalent 
to a point, and the Berry phase must be quantized as 
2irC, where C is the integer-valued Chern number. This 
is the topological invariant described above that underlies 
the remarkable quantization of the quantum Hall effect 
3, 18]. The Chern number has previously been measured 
through the conductance and chiral edge states of quan- 
tum Hall systems, and the total Berry phase associated 
with a Dirac point in graphene has been detected from 
transport [III [III and ARPES [2l| measurements. To our 
knowledge, the local Berry curvature has not been mea- 
sured directly, although its influence has been detected 
on a number of physical quantities (22[. It also plays 
an important role in the anomalous quantum Hall effect 
0, [H H and in semiclassical dynamics @, IM IH • 

Ultracold gases present an excellent opportunity to 
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study the Chern number and the topology of bands di- 
rectly. Optical lattices can impose periodic potentials, 
and ultracold gas experiments are highly controllable, 
tunable, and clean. Recent theoretical papers have shown 
how to detect the Chern number for certain tight-binding 
models in time-of- flight measurements [13, Sj- I n this 
paper, we propose a much more general method, based 
on semiclassical dynamics, that can be used to directly 
map out the Berry curvature over the Brillouin zone for 
any form of underlying lattice. 

In the absence of Berry curvature, a wave packet sub- 
jected to a constant external force will execute Bloch os- 
cillations. These have not been observed for bulk crys- 
talline electrons due to electronic scattering off lattice 
defects, but Bloch oscillations have been seen in other 
physical systems including semiconductor superlattices 
[29l | and ultracold gases [3(1 HH . 

In the presence of Berry curvature, the wave packet 
dynamics will be strongly modified. In principle, these 
effects may be measured directly from the real-space tra- 
jectory of such a wave packet, as has been previously 
proposed [32], [33|. However, in two dimensions, Bloch 
oscillations are complicated [33 - [36| . and it can be dif- 
ficult to disentangle the effects of Berry curvature from 
the usual effects arising from the group velocity. We pro- 
pose how this may be overcome using a "time-reversal" 
protocol. This will allow experiments to map the Berry 
curvature over the entire Brillouin zone. 

The structure of the paper is the following. First, we 
introduce the Berry curvature and describe its effects on 
the semiclassical dynamics. We then discuss the compli- 
cations of 2D Bloch oscillations and outline a protocol 
for how the Berry curvature may be mapped experimen- 
tally. We illustrate this discussion with numerical results 
for the Berry curvature and semiclassical dynamics for 
three interesting specific models: the asymmetric honey- 
comb lattice and two optical flux lattices [iHl [l7| . Finally 
we discuss general experimental considerations and time- 
of-flight experiments. 



II. THE BERRY CURVATURE AND 
SEMICLASSICAL DYNAMICS 

From Bloch's theorem, the eigenfunctions of a peri- 
odic potential can be expressed as tp n M{ r ) = Un,k(r)e lk ' r , 
where the Bloch function, Un^r), has the periodicity of 
the underlying lattice and n is the band index. We con- 
fine our discussion to the 2D xy plane, where the band 
structure is characterized by the energy, e n (k), and the 
scalar Berry curvature, fi„(k). All the concepts can be 
extended to three dimensions, where the Berry curvature 
must be treated vectorially. 

The Berry phase, "f n , for adiabatic transport in k space 
around a closed curve C bounding a region S is 



A„(k) = i(u nik \ — \u njk }, 
0„(k) = V k x A„(k) • z, 



(2) 
(3) 



where A„ (k) is the Berry connection, a gauge-dependent 
vector potential, and fl n {k) is the Berry curvature. The 
Berry curvature is a gauge-invariant, local property of 
the band structure. It vanishes when both time-reversal 
and inversion symmetries are present. The Berry phase 
is geometrical and similar to the Aharonov-Bohm phase, 
with the Berry curvature playing the role of a magnetic 
field. When integrated over the whole Brillouin zone, 
the resultant Berry phase is quantized and equal to 27r 
times the Chern number. The Chern number is thus a 
global topological property of the band, and underlies 
phenomena such as the integer quantum Hall effect |J. 
However, the local Berry curvature also has important 
physical consequences, for instance, on the semiclassical 
motion of a wave packet @, [22|, H^, [26| . 

To describe the semiclassical dynamics, we consider a 
gas of noninteracting fermions or bosons that is prepared 
in a wave packet with a center of mass at position r c 
and momentum k c |30L l37| . For atoms initially prepared 
in the bottom of the lowest band, the temperature con- 
tributes to the initial momentum spread of the atoms. 
We therefore assume that the temperature is less than 
the bandwidth so that the wave packet does not cover 
the whole Brillouin zone. We then consider a constant 
external force F, which in a solid state system would 
usually be due to an electric field. However, ultracold 
gases are neutral, and this force instead may come from 
linearly accelerating the lattice |30j, |38[ or from gravity 
[3ll Hs-HTJ. It is assumed that the force is sufficiently 
small that the motion is adiabatic such that Landau- 
Zener tunneling [42] is negligible and the wave packet re- 
mains in a single band; henceforward we drop the band 
subscript n. The semiclassical equations of motion are 
then HI 



I de(k c 
h dk r 



(kc x z)f2(k c ), 
hk c = F. 



(4) 
(5) 



7„ = f dk- A„ (k) = I dS- zfi„ 
Jc Js 



(k), (1) 



We shall refer to the first term in (j4|) as the group ve- 
locity and to the second as the Berry velocity. We note 
that these semiclassical equations are further modified if 
there is an external "magnetic" field [25] in addition to 
the external force F. We shall not discuss this further 
here, assuming that any magnetic field has the periodic- 
ity of the lattice and is incorporated into the bandstruc- 
ture through the (magnetic) Bloch states Uk( r ) (see Sec. 
IV Bjl . The effect of a magnetic field on Bloch oscillations 
has also been analyzed directly from the tight-binding 
Hamiltonian in Refs. |43| and 144 

To theoretically simulate the semiclassical dynamics, 
we must be able to calculate the Berry curvature. In 
general, a simple analytic expression is not possible and 
the Berry curvature is calculated numerically. This re- 
quires a discretized version of ©, as eigenfunctions are 
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found computationally over a grid in k space. There is 
an inherent phase ambiguity in the Bloch states, and so a 
gauge must be chosen to calculate the Berry connection. 
The Berry curvature is gauge invariant, and can be found 
on this grid by the method of Fukui et al. [45j , which ap- 
plies a geometrical formulation of topological charges in 
lattice gauge theory, where the Berry curvature is cal- 
culated from the winding of U(l) link variables around 
each plaquette in the Brillouin zone. We have used this 
method throughout this work for the numerical calcula- 
tions. 

It is of interest to note that effects of Berry curva- 
ture also arise in the semiclassical dynamics of a wave 
packet in a time-dependent one-dimensional optical lat- 
tice [22|, EH, |47| . The Berry curvature is then denned over 
a 2D parameter space made up of the one-dimensional 
quasimomentum and time. The Bloch oscillations of a 
wave packet in such a potential were theoretically stud- 
ied in Ref. |U 



The real-space Lissajous-like figures can be compli- 
cated two-dimensional oscillations, bounded by xbi oc 
VBiTBi (Fig. Q]). For the Berry curvature to change this 
trajectory significantly, it would be necessary to wait un- 
til the wave packet drifts outside of the bounding box. As 
a result, experiments would measure only the net Berry 
curvature encountered along a path. Information would 
be lost about how the Berry curvature is distributed in 
momentum space and notably whether its sign changes. 

Furthermore, in 2D there can be an additional drift 
in the wave packet's position, independent of the Berry 
curvature, if the wave packet does not start at high sym- 
metry points such as the zone center k = (0, 0) [3a.l4a|. 
Thus, merely observing a transverse drift is not, by itself, 
conclusive evidence of nonzero Berry curvature. 

IV. A "TIME-REVERSAL" PROTOCOL TO 
EXTRACT THE BERRY CURVATURE 



III. 2D BLOCH OSCILLATIONS 

The pioneering experiments on Bloch oscillations in ul- 
tracold gases were (quasi-)one dimensional [10, Hl|, and 
only recently has the extension to 2D been investigated 
10]. In 2D, Bloch oscillations have various interesting 
features in their own right, even before the Berry curva- 
ture is considered. 

One important consequence of dimensionality is that 
the real-space Bloch oscillations in 2D become Lissajous- 
like [3J-136J. For separable potentials, ID Bloch oscilla- 
tions along the x and y axes are simply superposed. For 
an arbitrary force F = [F x , F y ), the wave packet's motion 
is periodic along ki with periods TBi — h/\Fi\a (where i 
runs over x,y). The ensuing motion depends on the ratio 
F x : F y . For nonseparable potentials, studies show that 
similar behavior can be expected when the force an 
is weak and Landau-Zcncr tunneling is negligible [34 



)lied 
35]. 




FIG. 1. An example of a Lissajous-like figure for the square 
optical flux lattice (Sec. IVB1|) . The ratio F x : F y is 9 : 10, 
and a low force of |F| = 0.05Fr is used to minimize the ef- 
fects of Berry curvature over one oscillation (Fr = h 2 /2m\ 3 
where A is the optical wavelength and m is the mass of the 
atomic species; see Sec. [Vj. The Lissajous-like figure is ap- 
proximately bounded by the Bloch oscillation lengths, and so 
it obscures the effects of Berry curvature within this box. 



Berry curvature effects can be isolated by considering 
the dynamics under a reversal of "time." In doing so, it 
is important to be able to measure the velocity of the 
wave packet. We shall discuss in Sec. I VII how this may 
be done in experiments: for instance, through tracking 
the position of the wave packet in real space or through 
the momentum distribution, as in the seminal paper of 
Ben Dahan et al. [30j . 

We consider first measuring the velocity for a given 
force, F, at a particular point, k, in the Brillouin zone. 
This can be achieved in an experiment in which the wave 
packet has evolved according to k(t) = k(0)+Fi/7i. This 
velocity is uniquely defined (within the single band ap- 
proximation) at each point k along the trajectory and we 
denote it as Vk(+F). We now consider measuring the ve- 
locity in an experiment in which the wave packet passes 
through the same point k, but with opposite direction 
of the force, — F, which we denote Vk(— F). This can be 
achieved, for example, by evolving the wave packet along 
the line k(0) + Ft/Ti for a time T that moves it past the 
point of interest (e.g., to the Brillouin zone boundary), 
and then retracing this path using the reversed force — F. 
From ((3]), we can see that the Berry velocity changes sign, 
while the group velocity remains invariant. The two ef- 
fects can then be separated: 



v k (+F)-v k (-F) = --(FxzMk), 



v k (+F)+v k (-F) 



2 de(k) 
h <9k ' 



(0) 
(7) 



This transformation is equivalent to a time-reversal op- 
eration, and it cleanly removes the effects of the complex 
Lissajous-like figures in 2D. 

The Berry curvature can now be found at each point 
along the wave packet's path. By varying the path across 
the whole Brillouin zone, the Berry curvature is mapped 
out and the Chern number is directly measured. The 
path of the wave packet may be chosen in various ways, 
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FIG. 2. Two methods for mapping the Brillouin zone, (a) 
Rotating the force with respect to the lattice between exper- 
iments allows each wave packet to travel a different path, (b) 
With a large ratio F x : F y , a single wave packet successively 
travels many paths within the Brillouin zone. 



two of which are illustrated in Fig. [5J First, the align- 
ment of the force with the lattice may be rotated so that 
different trajectories are successively explored. In such a 
scheme, the measurement time can be short, correspond- 
ing to the time taken for the wave packet to travel once 
across the Brillouin zone. However, it would be impor- 
tant to align the force precisely each time. An alternative 
scheme would be to make the ratio F x : F y large. The 
wave packet will cover a large area of the Brillouin zone 
during a single Bloch oscillation. The force needs to be 
aligned only twice (for +F and — F) but longer measure- 
ment times would be required. A combination of these 
methods may be most suitable. 



A. Relation to the Chern number 

By this approach measurements of the velocity of the 
wave packet can be used to extract the Berry curvature. 
Measurements of the net drift of the wave packet in po- 
sition space can be used to measure the Chern number 
of the band. To illustrate the idea, it is convenient to 
consider a Brillouin Zone (BZ) of rectangular symmetry, 
and the set of paths k = (0, k y ) — > (K x ,k y ) that are 
traced out by a force +F in the x direction, for different 
values of the initial wave vector k y . Here K x denotes 
the size of the reciprocal lattice vector in the x direction, 
so the path traverses the full width of the Brillouin zone 
once. The net drift of the wave packet in the y direction, 
Ay = J Vydt, is 

r K * i de r Kx 

Ay{k y ,+F)= dk x + Q(k x ,k y )dk x . (8) 

Note that, for general k y , there is a transverse displace- 
ment not only from the Berry curvature but also from 
the group velocity [36|, IHj]. Reversing the force, such 
that the set of paths run in the opposite direction, over 
k = (0, k y ) — > (—if a; i k y ), the displacement becomes 



&y(k y ,-F) 



K 



* 1 de 
FdkZ 



K, 



Q(k x ,k y )dk x . (9) 



Thus, the contribution from the group velocity stays the 
same, but the contribution from the Berry curvature 



changes sign. Averaging the difference, 
Ay(k y ,+F)-Ay(k y ,+F) _ ' K 



Q{k x ,k y )dk x (10) 



extracts the part that depends on the Berry curvature. 
It is interesting to note that this contribution is inde- 
pendent of the magnitude of the force. The size of the 
transverse displacement is just set by the length scale of 
the underlying lattice (the lattice constant), and a nu- 
merical factor that involves the average Berry curvature 
along the trajectory. Furthermore, the total Chern num- 
ber can be found by summing over the set of trajectories 
with different values of k y , which can be used to represent 
an evaluation of the integral: 



C=±- 
2tt 



K y rK x 

dky I dk x £l{k x ,k y ) 
Jo 



(11) 



at discrete points in k y . (Clearly this approach can be 
readily adapted to a lattice of any symmetry, provided 
the set of paths spans the entire Brillouin zone once.) 

Note that, if instead of a wave packet, the band is filled 
(e.g., by noninteracting fermions) the Chern number may 
be measured from the net current density when a force, 
F, is applied: 



1 



Ky 



c 

dk x fl(k x ,ky)(F x z) = — (F x z). 

For a trapped gas, this result can be applied locally, with 
F set by the local potential gradient to give equilibrium 
currents. 



V. EXAMPLE SYSTEMS 

In this section, we illustrate our proposed method for 
measuring the Berry curvature for three example systems 
that are of experimental interest: the as ym metric hexag- 
onal lattice, and two optical flux lattices |16l Il7j for which 
the Chern number is nonzero. 

In optical lattices, the natural energy scale is set by 
the recoil energy, En = ft, 2 /2m A 2 , where A is the optical 
wavelength. Similarly, we can define a recoil velocity, 
v n = h/mX, and a unit of force, Fr = h 2 /2m\ 3 . Hence 
forward, we express all quantities in these units. 

The magnitude of the external force significantly af- 
fects the dynamics, as discussed further in Sec. IVI1 In 
previous experiments, the force has been introduced by 
linearly accelerating the lattice [3(J l38l. where the mag- 
nitude can be varied, or by gravity |3lL l39T - l4ll |. In our 
units, |mg| = 0.7-Fr for 174 Yb and for A = A « 578 
nm, the resonance wavelength coupling the ground and 
excited state in 174 Yb 14} . This choice of parameters is 



especially relevant to the optical flux lattices discussed 
below. We therefore primarily focus on the representa- 
tive case |F| = IF R . 
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For the evolution of the wave packet to be adiabatic, 
the rate of Landau-Zener tunneling to the next lowest 
band must be small. The probability of a Landau-Zener 
transition where the bands almost touch is given by [42j 



p = e- a " /a , 



(13) 



where a is the acceleration of the atoms moving under 
the external force, a c — (fe) 2 A/8?i 2 , and Ss is the band 
gap. This can therefore be neglected when the force is 
small or the band gap is large. 



A. The asymmetric hexagonal lattice 

The tight-binding hexagonal lattice has long been stud- 
ied in condensed matter physics as a simple model for 
graphene |49j. Thanks to recent advances, optical lat- 
tices of hexagonal symmetry (or closely related forms) 
can be imposed on ultracold gases and phenomena asso- 
ciated with the interesting band topology can be directly 
studied HSIM EH- 

In the presence of both inversion and time-reversal 
symmetry, the bands touch at two Dirac points in the 
corners of the hexagonal Brillouin zone. If either of these 
symmetries is broken, band gaps open and Berry curva- 
ture appears at these points, as in the Haldane model 
where time- reversal symmetry is broken [52j | . The Chern 
number has also been observed experimentally for time- 
reversal symmetry breaking in graphene [l9|, 0] • 

The asymmetric hexagonal lattice instead breaks in- 
version symmetry by introducing an onsite energy differ- 
ence between the two lattice sites, A and B. The opening 
of band gaps with asymmetry has already been studied 
experimentally in graphene [53j . but Berry curvature ef- 
fects have not been observed directly. Theoretically, the 
Berry curvature can lead to a quantum valley Hall effect, 
which may be useful for valley-based electronic applica- 
tions [22], [49|, HIJ ■ It would therefore be of great interest 
to study this system in ultracold gases. 




FIG. 3. The asymmetric hexagonal lattice in real space (on 
the left) and reciprocal space (on the right). The real space 
lattice vectors are ai = a(v3/2, — 1/2) and a 2 = a(0, 1), 
where a is the lattice spacing (and for this geometry, a — 
2A/3). The sublattices are connected by Ri = a(l/2v / 3, 1/2), 
R 2 = a(l/2vi},-l/2), and R 3 = a(-l/VE,0). The dotted 
lines indicate the unit cell. The reciprocal lattice vectors are 
then Ki = (4tt/a/3o)(1, 0) and K 2 = (4 7 r/ v / 3a)(l/2, V3/2). 



The honeycomb lattice can be viewed as two interpen- 
etrating triangular sublattices, for A and B, each with 
one site per unit cell (Fig. [3|). With on-site energies of 
±W on A/B sites, and including only nearest neighbor 
hoppings, the Hamiltonian is 



H(k) 



W V(k) 
V*(k) —W 



where V(k) = -J[e lkRl 
energy bands are then 



gik-R.2 _|_ gik-R.3 



(14) 



The two 



e(k) = ±y/W 2 + \V(k)\< 



(15) 



For W = 0, the energy bands have two Dirac 
points at which |V(k)| = 0: these are at k = 
(4tt/\/3o)(1/2,1/2VS) and k = (47r/V3o)(0, 1/V3), 
which we label as K and K' . Near each of the 
Dirac points, the effective Hamiltonian takes a sim- 
ple form. Close to the Dirac point K, writing k = 
(47r/v / 3a)(l/2, 1/2V3) + q, the effective Hamiltonian is: 



H(q) = 



W 

-hv F {q x ■ 



-hv F (q x - iq y ) 
-W 



(16) 



where Tivp = (V3/2)aJ: a Dirac equation with mass. 
The band structure is shown in Fig. Ufa) for W = 0.5Er 
and t — I.OEr. For this value of W, the band gap at the 
Dirac points is 1.0-Er and the Landau-Zener tunneling 
probability is less than 0.09 for |F| = IFr. Near the 
Dirac point, the Berry curvature is |22j | 



fi(q) 



2(W 2 



3/2 



(17) 



For the Dirac point K' , writing k = (47r/V^a)(0, 1/VS)+ 
q, the Berry curvature has the same form but opposite 
sign. 
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(b) 



-3 

k x a 



FIG. 4. (a) Band structure of the lowest band of the asym- 
metric hexagonal lattice for W — 0.5Er and t = I.OEr, and 
the energy in units of Eh (for this geometry, a = 2 A/3). Due 
to the asymmetry, gaps have opened at the Dirac points at 
the corners of the Brillouin zone, (b) The Berry curvature 
mapped out for W = 0.5Er, using the method of Ref. |45| . 
Light shading indicates Q, > and dark shading indicates 
Q < 0. Positive and negative regions cancel, giving a net 
Chern number of zero. 

The resulting map of Berry curvature for the asymmet- 
ric hexagonal lattice is displayed in Fig. HJb). This was 
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previously found analytically in Ref. [5ol The Berry cur- 
vature around points K and K' has opposite signs such 
that the net Chern number of the band is zero. This 
vanishing Chern number is required by the fact that the 
system is time- reversal invariant. 




FIG. 5. (a) The trajectory of a semiclassical wave packet 
through the Brillouin zone, starting from k = (0, 0) with F = 
(0.0, I.OFr). The numbers indicate the order in which the 
path is traveled, (b) The corresponding real-space trajectory 
of the wave packet, starting from the origin, for W = 0.5-En. 
This result was previously obtained in Ref. I33L The motion 
along x is due to the Berry curvature, while that along y is 
due to the band structure. 



that the 2D Bloch oscillation is also simply directed along 
y. This real-space trajectory was previously obtained in 
Ref. |H, where the effects of a perturbing "magnetic" field 
were also discussed. 

The velocities along x and y are shown in Fig. [6] As 
the wave packet passes through K', the negative Berry 
curvature gives it a positive velocity in the x direction. 
In between K' and K, there is no curvature and it moves 
with a group velocity along y. When it passes through 
K, the positive curvature gives it negative x velocity. As 
the regions of curvature have t he same magnitude, the 
effects cancel and there is no net drift. 

For F = (0, F) it is simple to determine the Berry cur- 
vature because the group and Berry velocities are per- 
pendicular. For more general directions of the force, 
Lissajous-like oscillations will make it difficult to extract 
any information about the Berry curvature from the real- 
space motion. 

As proposed above, the Berry curvature may be cleanly 
mapped from the velocities using a time-reversal proto- 
col. This is illustrated in Figs. [5tc) and^d). Here the 
velocities along x and y for +F and — F are combined to 
show the Berry velocity and the group velocity, respec- 
tively. 




-4 4 .4 4 




-4 4 -4 4 



FIG. 6. (a) Velocity along a; of a wave packet moving under 
the force F = ±(0.0, I.OFr) with W = 0.5E R . k y is measured 
along the path travelled (which is periodic in 47r). (b) Velocity 
along y. (c) Applying the time-reversal protocol to extract 
the Berry velocity: v(+F) - v(-F) = -2/h(F x z)fi(k). 
(d) Applying the time-reversal protocol to extract the group 
velocity: v(+F) + v(-F) = (2/fi)de(k)/dk 



B. Optical flux lattices 

One of our main motivations for mapping the Berry 
curvature is to find a way of experimentally characteriz- 
ing optical flux lattices. These have lately been proposed 
as schemes to access fractional quantum Hall physics in 
ultracold gases 0, [l7|. The optical flux lattices have 
bands with nonzero Chern numbers. They consist of a 
state dependent potential in register with an interspecies 
coupling, and lead to effective magnetic flux with a high 
nonzero average (with a nonzero integer number of mag- 
netic flux quanta per unit cell). The specific optical 
coupling depends on the geometry and atomic species 
chosen. For atoms with a ground state and long-lived 
metastable excited state, such as an alkaline earth atom 
or ytterbium, a s imp le one-photon coupling scheme can 
be implemented [161 ]. For more commonly used atomic 
species, such as alkali atoms, two hyperfine states can be 
used with coupling via two photon processes [17| . Here, 
we discuss an example of each: the square one-photon 
optical flux lattice and the F = 1/2 two-photon optical 
flux lattice. The lowest energy band in both lattices can 
have a Chern number of one and hence is topologically 
equivalent to the lowest Landau level, allowing quantum 
Hall physics to be accessed. 



From the Berry curvature and band structure, we can 
now calculate the semiclassical motion of a wave packet 
in this system (Fig. [3]). To illustrate clearly the effects of 
Berry curvature, we start the wave packet at k = (0, 0) 
and consider a force aligned along the y direction, so 



1, One-photon square optical flux lattice 

In this scheme, the electronic ground state and a 
long-lived excited state are coupled via a single photon 
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process [14|. The Hamiltonian in the rotating wave ap- 
proximation is 



„ = S-i + V< r) , 



(18) 



where 1 and V are 2x2 matrices acting on the two 
internal states of the atom. We neglect interactions, an 
assumption that is discussed further in Sec . I VII The 
square optical flux lattice is generated when [16| 

V sq = V[cos(nx)a x + cos(ny)a y + sin(nx) sm(Ky)a z ], 

(19) 

where V sets the energy scale of the potential, <f,; are 
the Pauli matrices, k = 2n/a, and the lattice vectors are 
ai = (a, 0), &2 = (0, a). The flux density is everywhere 
of the same sign, and leads to a total flux per unit cell of 
= 2 El. 
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FIG. 7. Band structure of the lowest band in the magnetic 
Brillouin zone for the one-photon square optical flux lattice 
with (a) V = 0.0E R , (b) V = 1.0E R and (c) V = 3.0E R , 
and the energy in units of Er. At V = O.OEr, the lowest 
band touches the second lowest band at k = (n/a,±ir/a). As 
V is turned on, band gaps open at these points and Berry 
curvature is formed. 



which do not commute but satisfy T 2 T\ — —T\Ti- These 
operators represent rotations in spin space and transla- 
tions by jai 2, which enclose half a flux quantum (as 
N<f, — 2) [TfJ,[l8(. The magnetic Brillouin zone is defined 
by a unit cell containing an integer number of flux |18| ; we 
choose a cell containing a single flux with vectors ai and 
a2/2. The corresponding commuting operators are 
and T 2 , with eigenvalues e lk ai and e lk a2 / 2 . This defines 
the Bloch wavevector, k, and the associated magnetic 
Brillouin zone Hi . Now the first Brillouin zone covers 
—ir/a < k x < ir/a and —2i:/a < k y < 2n/a, doubling in 
size. Thanks to this unfolding, the lowest band is nondc- 
generate at each Bloch wavevector k. The resulting band 
structure is shown in Fig. [7] for V = 0.0E R , V = 1.0E R , 
and V = 3.0E R . 

The Berry curvature is shown over the magnetic Bril- 
louin zone for the lowest band in Fig. [8l for V — I.OEr 
and V = 3.0Er. For nonzero V, the Chern number of 
this band is one, so it is analogous to the lowest Lan- 
dau level. For small V the Berry curvature, Q, is highly 
peaked at positions k = (n/a, ±7r/a); as V increases, f2 
spreads out while remaining centered on these two points. 






36.7 
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0.0 





FIG. 8. Contour maps of the Berry curvature, f2(k), for 
the square optical flux lattice when (a) V = I.OEr and (b) 
V = 3.0E R . n is in units of 1/(MBZ area) = a 2 /8vr 2 . SI is 
concentrated at points k — (n/a, ±7r/o), spreading out with 
increasing V . 



We consider the regime where V < h 2 K 2 /2m = AEr 
(for this geometry a = A/2). We expand the periodic 
Bloch functions over a set of reciprocal lattice vectors K, 

1 / l(n,k) \ 

and diagonalize the resulting Hamiltonian to find the 
band structure. As we consider low V, a small, finite 
set of K vectors will give the eigenfunctions and values 
to within the required numerical accuracy. 

The eigenfunctions are everywhere twofold degenerate, 
corresponding to two magnetic subbands. To distinguish 
between these states, we reinterpret the system within 
the magnetic Brillouin zone (MBZ) [II [lj, [I! . The op- 
tical coupling is invariant under the magnetic translation 
operators: 

f! = V 4ai,V f 2 = a x e^ v (21) 



In Fig. [Hlb) we illustrate the real-space trajectory for 
semiclassical motion of a wave packet in this optical flux 
lattice, for V = OAEr and V = 3.0Er. As discussed 
above, we take |F| = I.OFr, which is representative of 
the gravitational force |F| = rag for 174 Yb and A = 578 
nm, parameters which are particularly relevant for this 
scheme [3 [3. As shown in Fig. in momentum 

space the wave packet starts at k = (0,0) and moves 
under a force F parallel to the (1,1) direction such that 
it passes through the points k = (7r/o, ±7r/a) at which 
there is large positive Berry curvature. With F aligned 
along (1,1), the group velocity is parallel to the force, 
while the Berry velocity is perpendicular (Fig. [TU)) . As 
the Berry curvature is everywhere positive, there is a net 
drift as successive regions of high SI are traversed. Fol- 
low potentials such as V — OAEr and V = 3.0Er, the 
bandgap is small and the probability of Landau-Zener 
transitions on crossing the Brillouin zone boundaries is 
large. In an experiment where |F| = rag, this probabil- 
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ity can be reduced below 0.1 by increasing the potential 
above V — 3.2E R . With higher V, the Berry curvature 
spreads out and hence the trajectory bends along more 
of its length. For clarity, we therefore discuss V = OAEr 
where the effects of Berry curvature and the group ve- 
locity are easiest to understand. The dynamics will be 
qualitatively the same for higher V, and we discuss the 
dependence of the motion on both V and |F| in more 
detail in Sec. [HI 
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FIG. 9. (a) The trajectory of a semiclassical wave packet 
through the magnetic Brillouin zone, starting from k = (0, 0) 
with |F| = l.OFfl directed along (1,1). The numbers indicate 
the order in which the path is traveled. The dotted line indi- 
cates the simple Brillouin zone, while the dashed lines shows 
the extension into the magnetic Brillouin zone, (b) The cor- 
responding real-space trajectory of the wave packet starting 
from the origin, for V = OAEr and V = 3.0Fr. The motion 
perpendicular to F is due to the Berry curvature, while that 
parallel is due to the bandstructure. 
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FIG. 10. The velocity of the wave packet (a) perpendicular 
and (b) parallel to the applied force for |F| = I.OFr directed 
along (1,1) and V = OAEr. k is measured along the path 
travelled. For this simple alignment of the force, (a) contains 
the effects of Berry curvature along the path, while (b) shows 
the group velocity. 

As before, the simplicity of the trajectories and veloci- 
ties relies on the alignment of F along a special direction. 
More generally in 2D, the Bloch velocity will not be par- 
allel to the force and complex Lissajous-like figures will 
be observed; an example of this was previously shown in 
Fig. [T]for one oscillation at a low force, |F| = 0.05Fr, 
where the effects of Berry curvature are small. 

Here we demonstrate the effect of Lissajous-like figures 
on the motion for parameter ranges of interest. Figure [TT1 



shows the real-space trajectories for (a) |F| = I.OFr and 
(b) when the force is aligned such that F x : F y — 

1 : 16. The Bloch motion is no longer purely parallel to 
the force and while an average drift may be measured, 
the details of the motion due to Berry curvature have 
been lost. 
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FIG. 11. Trajectory for a wave packet traveling with (a) |F| = 
I.OFr and (b) |F| = 5.0Fr where the force is aligned such 
that F x : F y = 1 : 16 and V = 0AE R . The Bloch motion 
is no longer purely along the direction of the force, and the 
resultant trajectory is complex. The net drift between Bloch 
oscillations is a measure of the total Berry curvature along a 
path, but other information is obscured. The trajectory due 
to the group velocity increases in size relative to the trajectory 
due to the Berry curvature with decreasing force (Sec. IVI|) . 
leading to more complicated motion in (a) than (b) . 
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FIG. 12. (a) Velocity along a; of a wave packet moving under 
the force |F| = I.OFr with F x : F y = 1 : 16 and V = 0AE R . 
k is measured along the path traveled. (b)Velocity along y. 
(c) Applying the time-reversal protocol to extract the Berry 
velocity. The Berry velocity calculated from v x and v y differ 
by a factor of 16 from the ratio of the forces, and by a sign, 
due to the cross product, (d) Applying the time-reversal pro- 
tocol to extract the group velocity. As expected, there are 16 
oscillations from v y for every one from v x . 



We illustrate how our time-reversal protocol may be 
used to extract the local Berry curvature. Figures Q2] 
and [T3] show the velocities along x and y over the path 
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FIG. 13. As in Fig. [IS] but now for |F| = 5.0F R . 



through the Brillouin zone. By comparing the time- 
reversed velocities, quantities proportional to the Berry 
velocity and the group velocity are extracted. Figure [T4l 
shows the resulting map of Berry curvature over the Bril- 
louin zone for |F| = 1.0-Fr. The same result is obtained 
(up to a scale factor) from either of these two cases of 
different magnitude of the force. 




FIG. 14. \v x (+F) - v x (—F)\/2 = l/h(F y n) plotted along 
the path taken by the wave packet in the Brillouin zone for 



|F| = 1.0F H with F x 



1 : 16 and V = 0AE R . This 



technique will enable experiments to directly map out the 
Berry curvature. 



2. Two-photon optical flux lattice for F = 1/2 

To generate optical flux lattices for the more commonly 
used alkali atoms one must employ dressed states involv- 
ing a two-photon coupling [17]. We consider the rep- 
resentative case of a lattice with triangular symmetry 
where the hyperhne states coupled have angular momen- 
tum F = 1/2, as for 171 Yb. Qualitatively similar semi- 
classical dynamics are expected for other values of F, 



such as the experimentally common case of F = 1 for 
87 Rb [17|. The two-photon optical flux lattice we study 
leads to a net effective magnetic field in real space in 
which there is one flux quantum per unit cell. Semiclas- 
sical motion within a similar scheme has previously been 
studied in Ref. HH In that scheme, the artificial mag- 
netic field is still locally nonzero, but the flux per unit 
cell vanishes [171 |32| and the Chern number of each band 
is zero. For the optical flux lattice described here, the 
energy bands may have nonzero Chern numbers. 

For the F = 1/2 optical flux lattice, two hypcrfinc 
ground states, g±, with angular momentum J g — 1/2, are 
coupled to an excited state, e, also with angular momen- 
tum J e = 1/2, via an off-resonance excitation that en- 
sures the population of e remains negligible. The Hamil- 
tonian then acts in the g± manifold, with the form of 
(IT8l) . The details of the optical cou pling and the resulting 
Hamiltonian are discussed in Ref. [171 ] and in Appendix 

m 

The final Hamiltonian can be written as: 

H> = &HU = (p ~^ k3/2)2 + V', (22) 
2m 

where k3 = fc(0, 1,0), and U is a unitary transforma- 
tion that is applied to expose the full symmetry of the 
system. The transformed optical potential V' has the 
maximal translational symmetry, causing p/h to be con- 
served up to the addition of the reciprocal lattice vec- 
tors Ki/2 = -fc/2(±V3,3,0) (Appendix©. The result- 
ing Brillouin zone, defined by K^ 2 , is equivalent to the 
asymmetric hexagonal lattice in Sec. IV Al As can be 
seen, an important feature of the two-photon optical flux 
lattice is that, under this unitary transformation, the mo- 
menta of g± are offset by ±^3/2. This offset does not 
affect the semiclassical equations of motion, which de- 
termine the rate of change of (crystal) momentum under 
an applied force ((5|), and which still apply for the bands 
formed from the eigenstates of H' . 

The optical coupling, V, is characterized by parame- 
ters e and 9 as well as the overall strength of the potential 
V (Appendix |A"]) . When 9 = e = 0, the optical potential 
does not couple the states g+ and g_, and acts on each 
simply as a scalar potential, with the same symmetries 
as the hexagonal lattice discussed above when inversion 
symmetry is unbroken. The offsets of the momenta of g± 
shift the Dirac points of these two states relative to each 
other. A small nonzero 9 breaks inversion symmetry and 
opens up gaps at the Dirac points in such a way that the 
pairs of bands are topologically trivial, i.e., have a net 
Chern number of zero. Nonzero e and 9 together break 
time-reversal symmetry and lead to bands with nonzero 
Chern number. 

We have numerically calculated the band structure and 
Berry curvature of the lowest band for a case of large cou- 
pling, when V = 1.8E R , 9 = 0.3, and e = 0.4, for which 
this band has a Chern number of one (Fig. [T5|) . Note 
that in this case, the energy minimum of the band is not 
at k = (0,0). Applying the force along (0,1), the wave 
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FIG. 15. (a) The energy (measured in Er) of the lowest band 
for V — 1.8Er, 9 — 0.3, and e = 0.4. For these parameters, 
the lowest band has a Chern number of one. (b) The corre- 
sponding Berry curvature. fi(k) is in units of 1/(BZ area). In 
this case, the Berry curvature is significant over much of the 
Brillouin zone and has regions of both positive and negative 
sign. 




x / a 

FIG. 16. The real space trajectory of a wave packet under a 
force F = (0.0, O.LFr) over one period, for V = 1.8E R , 9 = 0.3 
and e = 0.4. The wave packet follows the path described in 
Fig. [5] While the group velocity is again purely along y, the 
motion due to the Berry curvature is more complicated due 
to the varying sign of Q.. 

packet will follow the same path as in the asymmetric 
hexagonal lattice (Fig. [5]). We choose |F| = Q.IFr to en- 
sure that the Landau-Zener tunneling probability is small 
(below 0.1). The resulting real-space trajectory is shown 
in Fig. 1161 Due to the simple alignment of the force, we 
can again associate the y motion with the group velocity 
and the x motion with the effects of the Berry curvature. 

As the Berry curvature is substantially spread out 
through the Brillouin zone (Fig. [T5|i . there is now a con- 
tinual drift along x over the trajectory. Since the Chern 
number is one, the Berry curvature is (largely) of the 
same sign along the trajectory, leading to a net drift of 
the wave packet over one period of the Bloch oscillation. 
In this lattice, there are also regions of both positive and 
negative Berry curvature. When the Berry curvature is 
positive, the wave packet travels in the negative x direc- 
tion, while when f2 is negative, the wave packet moves 
along the positive x axis. The sign change will therefore 
not be detected if only the total drift is measured. In- 
stead, the time-reversal procedure described above can 
again be applied to cleanly map out the local Berry cur- 
vature. 



VI. EXPERIMENTAL CONSIDERATIONS 

We now consider how the velocity may be measured 
and how feasible it will be to observe the Berry curvature 



effects experimentally. 

The mean atomic velocity may be extracted from the 
time-of-flight expansion image. This measures the mo- 
mentum distribution [56| . from which the mean mo- 
mentum (p) may be deduced by the weighted average. 
The mean atomic velocity of the initial wave packet 
then follows from Ehrenfest's theorem as (v) = (p)/m. 
This approach was successfully experimentally applied by 
Ben Dahan et al. [30( to detect Bloch oscillations in a 
one-dimensional lattice. The same approach applies for 
dressed states of internal atomic states, governed by (fl"51) ; 
in this case the mean velocity can be obtained from the 
average momentum over all internal states. 

Alternatively, it is possible to extract the velocity di- 
rectly from measurements of the center of mass motion 
in real space. Thanks to recent experimental advances, 
the position of the wave packet's center of mass may be 
imaged with a high resolution, on the order of a lattice 
spacing (57j . For bands with nonzero Chern number, the 
Berry curvature can cause the wave packet to have a net 
drift in space over each period of the Bloch oscillation; 
this leads to large cumulative effects on the position of 
the wave packet over many oscillation periods. As de- 
scribed in Sec. IIV Al measurements of the position of 
the wave packet therefore offer a sensitive way to show 
that the Chern number of the band is nonzero. Indeed 
measurements of these drifting trajectories in constant 
applied force are equivalent to measurements of the edge 
states that must arise for a band with nonzero Chern 
number. 

To this end, we consider how to maximize the im- 
portance of the effects of the Berry curvature relative 
to those of the group velocity. In the cases described 
above where the bands almost touch, we can consider 
the band gap, Se, as a small parameter. This applies to 
the hexagonal lattice when asymmetry is small and to the 
square one-photon optical flux lattice for small potential, 
V. (Note that the two-photon optical flux lattice is con- 
sidered far from the band closing regime.) In the small 
band gap limit, the momentum width 5k, over which the 
band is changed will be 5k ~ 5s/(hi)ji), where vr is the 
typical group velocity at the zone boundary for vanish- 
ing band gap. The Berry curvature is therefore nonzero 
over the area, A ~ (5k) 2 . Assuming the Berry curva- 
ture is uniform within this, the invariance of the Chern 
number means f2 ~ 1/A ~ l/(5k) 2 . The Berry velocity is 
vq = (k x z)f2 so, as the wave packet traverses this region 
in one Bloch oscillation, the Berry curvature leads to a 
displacement of size J v^dt ~ 5kQ ~ l/Sk ~ hvu/5e. 
Conversely, the typical group velocity is ~ vr so over the 
period of one Bloch oscillation, tb — h/(Fa), the typical 
amplitude of displacement is vrtb — hv^J (Fa). Thus, 
these contributions to the real-space trajectory have dif- 
ferent dependences on 5s: the Berry curvature contri- 
bution scales as l/5s, while the contribution from the 
group velocity is independent of 5s. This is found in our 
numerical results, but is shown only qualitatively in the 
results presented in Fig. \§\ as in this case 5s is not that 
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small. The effects of Berry curvature can therefore be 
maximized with respect to those of the group velocity by 
choosing a small band gap. Note also that the two con- 
tributions have different dependences on the size of the 
force: the displacement due to the Berry curvature is in- 
dependent of force, while that due to the group velocity 
is inversely proportional to it. Therefore, the effects of 
Berry curvature will be most evident for a high external 
force. Note, however, that there are some practical lim- 
itations on both the choice of force and band gap. To 
ensure that the evolution of the wave packet is adiabatic, 
the rate of Landau-Zener tunneling to the next lowest 
band should be negligible [Eq. [13] • The assumption of 
adiabatic evolution is therefore violated when the applied 
force is too high and the band gap is too small. Also, for 
small 5s, when the bands nearly touch, the Berry curva- 
ture becomes concentrated in small regions, of area (5k) 2 . 
This then requires the momentum of the wave packet 
and the alignment of the force to be precisely controlled 
in order to direct the wave packet through this region. 
For intermediate band gaps, the curvature is spread out. 
A natural compromise is to choose the band gap such 
that 5k is as small as the momentum uncertainty 1/w 
with which a wave packet can be prepared (here w is the 
spatial width of the wave packet); one then expects the 
displacement due to the Berry curvature to be l/5k ~ w 
on the order of the spatial extent of the wave packet. 

From our numerical calculations, we can quantitatively 
estimate the lengthscales of the dynamics. For example, 
we consider the dynamics of 174 Yb atoms in the one- 
photon optical flux lattice, with F = mg along (1,1), A = 
578 nm, and V = 3.2Er. For this choice of parameters, 
the Landau-Zener tunneling probability given by (|13l) is 
approximately 0.1. The wave packet follows the Brillouin 
zone path in Fig. [9ja) and has a real-space trajectory 
similar to that of the dashed line in Fig. HJb). As the 
wave packet moves along section 1 of its path [from k = 
(0, 0) to k = (ir/a, 7r/a)], it moves with an average group 
velocity of 1.0 mm s , traveling approximately 0.6 /zm 
in the direction of the force in real space. For V = 3.2Er, 
the Berry curvature is substantially spread out over the 
Brillouin zone, and so the wave packet's trajectory bends 
as it travels, moving it 0.3 /jm perpendicular to the force. 
As the wave packet continues from k = (— 7r/a,7r/a) to 
k = (0, 2tt/o), the group velocity changes sign and the 
wave packet travels 0.6 /im in the opposite direction to 
the force. The Berry velocity does not change sign, and 
so the wave packet moves a further 0.3 /iin perpendicular 
to the force. This behavior repeats for sections 3 and 4 of 
its path. The average Berry velocity over one complete 
oscillation is therefore approximately 0.3 mm s _1 . If the 
force is slightly misaligned, the trajectory will instead be 
a Lissajous-like oscillation, approximately bounded by 
a box of diagonal length 0.6 /im. For the same Berry 
velocity, the wave packet would then take approximately 
2.0 ms to drift this distance. These length and time scales 
are within current experimental capabilities. 

Two additional practical considerations are the effects 
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FIG. 17. (a) Oscillation in expectation value of a x for the one- 
photon optical flux lattice, as the wave packet travels along 
(1, 1). The dashed line is for V = O.OAEr, while the smooth 
line is V = OAEr. The area over which the transfer occurs 
decreases with potential in a similar way to how the Berry 
curvature area decreases, k is measured along the path of the 
wave packet, (b) Oscillation in population in the two states 
for the two-photon optical flux lattice as the wave packet trav- 
els along (0, 1) for V = 1.8E R , 6 = 0.3, and e = 0.4. The 
maximum population imbalance appears to be around points 
K and K'. k y is measured along the path of the wave packet. 



of dispersion and interactions. The wave packet will 
spread as it travels, and this dispersion could obscure 
the dynamics described. However, provided the center 
of mass of the wave packet can be measured to an accu- 
racy greater than its width, this should not prevent the 
observation of Berry curvature effects. 

Interactions destroy the coherence of a wave packet 
over time and can have a strong dephasing effect on Bloch 
oscillations [3l], US H3- Nonlinearity can also lead to 
the collapse of the wave packet into discrete solitons [HI, 
H§|. We have ignored the effects of interactions in our 
analysis, an approximation suitable over these time scales 
for fermionic atoms (39|, for species with low scattering 
lengths [4l]], or where the interaction strength can be 
tuned to zero by means of a Feshbach resonance [10, [6(| ■ 

One can also look for distinct features in the momen- 
tum distributions of atoms undergoing Bloch oscillations 
in the optical flux lattices discussed above (Fig. [T7"|) . 
In the two-photon optical flux lattice, the population of 
atoms oscillates between the two internal states as the 
wave packet moves through the Brillouin zone. Near the 
points K and K' , the population imbalance is maximum 
but of opposite sign, reflecting the characteristics of the 
Berry curvature. In the one-photon optical flux lattice, 
the unfolding into the full magnetic Brillouin zone means 
the Bloch states are eigenstates of & x . As a result, the os- 
cillation takes place between those superpositions of the 
internal states that are eigenstates of a x . The transfer 
between eigenstates occurs over an area which decreases 
with potential, in a similar way to how the Berry curva- 
ture area decreases. 



VII. CONCLUSIONS 

In summary, we have proposed a general method for 
mapping the local Berry curvature over the Brillouin zone 
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in ultracold gas experiments. The Berry curvature cru- 
cially modifies the semiclassical dynamics and so affects 
the Bloch oscillations of a wave packet under a constant 
external force. In particular, the Berry curvature may 
lead to a net drift of the wave packet with time. How- 
ever, two-dimensional Bloch oscillations are interesting 
in their own right, and one may lose information about 
the Berry curvature due to the complicated Lissajous-like 
figures that may arise. 

We have shown that this information can be recovered 
via a time-reversal protocol. The group velocity at a 
point in the Brillouin zone is invariant under a reversal 
of force, while the Berry velocity changes sign. As a 
result, the velocities under positive and negative force 
can be compared to extract either one or the other. This 
protocol will allow the local Berry curvature to be cleanly 
mapped out over the path of the wave packet through the 
Brillouin zone. 

We have demonstrated this protocol using the semi- 
classical dynamics of three model systems which are 
currently of experimental interest: the asymmetric 
honeycomb lattice and two optical flux lattices. Finally, 
we have discussed various experimental considerations, 
including how the velocity may be measured and how 
to maximize the magnitude of the Berry curvature 
effects on the dynamics. These methods can be used to 
characterize the topological character of band structures 
of complex optical lattices including optical flux lattices. 
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Appendix A: Hamiltonian for the Two-Photon 
Optical Flux Lattice for F — 1/2 

The Hamiltonian acts in the g± manifold, with the 
form of (fT5)l . The optical potential, V"(r), is characterized 
by the Rabi frequencies K m , m = 0, ±1, where mh is the 
angular momentum along z gained by the atom when it 



absorbs a photon (l7| . 

The potential is formed from laser beams at two fre- 
quencies: ujl and ujl + 5 where 6 is the Zeeman splitting 
of the ground states. The laser beams at wl are linearly 
polarised and traveling in the xy plane, while the laser at 
w l + 8 gives a plane wave propagating along the z axis 
with a z polarization. The resulting potential takes the 
form: 

,y_ hK tott, & (\k-\*-\k+\ 2 Ekv \ 

V ~ 3A i+ 3A ^ E K * |«+| 2 -|«-rV 

(Al) 

where k\ o1 = J2 m \ K m\ 2 , A = ujl — Ua, with to a as the 
atomic resonance frequency and it is assumed that | A | 3> 
\S\, \K m \. The field of the laser at frequency ujl + S is 
characterized by E, which serves as a uniform, adjustable 
coupling. 

The laser field at frequency ujl is formed from the su- 
perposition of three traveling plane waves of equal in- 
tensity and wave vectors in the xy plane. The set- 
up discussed in Ref. [17| has triangular symmetry, with 
an angle of 27r/3 between the beams. The wave vectors 
arc then k x = -fc/2( % /3, 1, 0), k 2 = fc/ 2 (V3,-l,0) and 
k3 = fc(0, 1, 0). Up to a scale factor and rotation, this is 
the same geometry as the asymmetric hexagonal lattice 
in Sec. IVA"1 

Each beam is linearly polarized at an angle 6 to the z 
axis, giving 

3 

k = k VV k< ' r [cos0z + sin0(z x kj)], (A2) 

where k is the Rabi frequency of a single beam. The 
relative strength of the laser fields at frequencies ui and 
+ S will henceforth be denoted as e = E/k. 
To define the unit cell, we consider a unitary transfor- 
mation: U = exp(— ik3 • rer 2 /2). This gauge transforma- 
tion is state dependent and leads to the Hamiltonian in 
Eq. ||22J), where V' = U^VU has the same form as V, 
with Ko replaced by k' — e~ lk3 ' r Ko. The coupling then 
only includes the momentum transfers Kj. 2 = ki 2 — k.3. 
These define the reciprocal lattice vectors of the largest 
possible Brillouin zone, and the smallest possible real- 
space unit cell. This unit cell is that of the hexagonal 
lattice discussed in Sec. IV Al 
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